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Computer tomography method for periodically moving objects 



The invention relates to a computer tomography method, in which a 
periodically moving object, especially a heart, is irradiated by a beam bundle. The invention 
relates in addition to a computer tomograph for carrying out the method and to a computer 
program for controlling the computer tomograph. 

In known methods of the kind mentioned at the beginning, the three- 
dimensional spatial course of the absorption or the attenuation of the radiation in the 
periodically moving object is reconstructed from measured values that are acquired using a 
detector unit. During that process, the effect of the periodic movement of the object is that the 
measured values contain information from different object states, which in turn leads to 
movement artifacts in the reconstructed data record. 

In order to reduce these movement artifacts, during the acquisition of 
measured values in known methods, movement signals, for example an electrocardiogram, 
dependent on the movement of the object, are recorded, so that a movement signal can be 
allocated to each measured value. Only measured values to which identical or only slightly 
differing movement signals are allocated are then used for the reconstruction. 

The disadvantage of this method is that identical movement signals do not as a 
rule correspond to identical object states. The measured values used for reconstruction may 
have been allocated substantially identical movement signals, but this does not mean that 
they have also been allocated substantially identical object states. The above-described 
selection of the measured values used for reconstruction therefore also leads to very marked 
movement artifacts. 



It is an object of the present invention to specify a computer tomography 
method, a computer tomograph and a computer program, with which these movement 
artifacts are less marked. 
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That object is achieved in accordance with the invention by a computer 
tomography method having the following steps: 

a) generation by a beam source of a beam bundle passing through a periodically 

moving object, 

5 b) generation of a relative movement between the beam source on the one hand 

and the object on the other hand, which comprises a rotation about an axis of rotation, 

c) acquisition by means of a detector unit, during the relative movement, of 
measured values that are dependent on the intensity in the beam bundle on the other side of 
the object, an acquisition time being allocated to each measured value and to the beam 

1 0 causing the respective measured value, 

d) detection of a movement signal depending on the movement of the object by 
means of a movement-detection device and determination of periods of the periodic 
movement by means of the detected movement signal, 

e) reconstruction of a computer tomography image of the object from the 

15 measured values, wherein there are used only measured values whose acquisition times lie 
within the periods (Ti...T 7 ) at time intervals (Ati... At 7 ; At'i... At' 7 ), which are so 
determined that a similarity measure applied to intermediate images of the same subregion 
(23; 25) of the object is minimized, wherein different intermediate images are reconstructed 
using measured values from time intervals (Ati . . . At 7 ; At* i . . . At' 7 ) from different periods 

20 (T|...T 7 ). 

Intermediate images of one and the same subregion are thus reconstructed 
using measured values from time intervals from different periods. That is, in each case 
exactly one period can be allocated to each intermediate image. The time intervals in the 
individual periods are now adjusted in such a way that, after a reconstruction of the 

25 intermediate images with the adjusted time intervals, a similarity measure applied to the 

intermediate images of the same subregion is minimized. This method can be applied to one, 
several or all subregions of the object that are reconstructable using measured values from 
time intervals from different periods. Finally, a computer tomography image is reconstructed, 
wherein exclusively measured values from the adjusted time intervals are used. 

30 The movement artifacts described in the beginning occur particularly 

noticeably in the regions of the reconstructed image that were reconstructed using measured 
values from time intervals from different periods. That is why it is important especially in 
these regions to use for the reconstruction measured values that were acquired during the 
time the object was in object states as identical as possible. This is achieved in accordance 
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with the invention by determining the corresponding time intervals so that a repeated 
reconstruction of the relevant subregion, in each case using measured values from different 
periods, leads to images as similar as possible. The subsequent reconstruction of the 
computer tomography image using all measured values that lie in the time intervals thus 
5 determined leads to a computer tomography image (CT-image) having fewer marked 
movement artifacts compared with the prior art. 

The term "periodic movement" is not restricted to a periodicity in the exact 
sense, that is, it is not restricted to movements in which object states regularly recur exactly, 
i.e. exactly identical object states at exactly equidistant instants. A periodic movement within 

10 the scope of the invention includes in particular movements that differ from a mathematical 
exactness, as is known in the case of periodically moving body organs, for example, the 
heart. That is, similar, substantially identical object states are passed through at instants that 
as a rule are substantially equidistant. 

Subregion within the meaning of claim 1 can also be a region to be 

15 reconstructed, which comprises the entire object or just a portion of the object. 

Claims 2 and 3 described methods in which each two consecutive periods 
form a period pair and a first intermediate image of a subregion is reconstructed exclusively 
using measured values whose acquisition instants lie in the time interval of the earlier period 
of a period pair, and a second intermediate image of the same subregion is reconstructed 

20 exclusively using measured values whose acquisition instants lie in the time interval of the 
later period. The intermediate images are compared with the similarity measure. The 
reconstruction of one of the intermediate images is now carried out with modified, especially 
shifted, time intervals until a break-off criterion dependent on the similarity measure is 
satisfied. This leads to a further reduction in movement artifacts. 

25 Claims 4 and 5 describe a break-off criterion and a similarity measure, which 

produce a reconstructed CT image of very good quality. 

In the embodiment as claimed in claim 6, each measured value is weighted 
before the reconstruction, which leads to a further improvement in the image quality. 

Claims 7 and 8 describe computer tomography methods in which use is made 

30 of a filtered back-projection respectively a lower three-dimensional resolution in the 

reconstruction of the intermediate images than in the reconstruction of the CT image, thus 
reducing the computing expenditure. 
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Claims 9 and 10 describe embodiments of the method according to the 
invention that use an electrocardiogram as movement signal. This leads to a CT image of a 
beating heart with only few movement artifacts. 

A computer tomograph for carrying out the method according to the invention 
5 is described in claim 1 1 . 

Claim 12 defines a computer program for controlling a computer tomograph as 
claimed in claim 11. 



10 These and other aspects of the invention are apparent from and will be 

elucidated with reference to the embodiments described hereinafter. 
In the drawings: 

Fig. 1 shows a computer tomograph, with which the method according to the 
invention can be implemented, 
15 Fig. 2 is a flow chart of the method according to the invention, 

Fig. 3 shows a schematic representation of a correlation between periods, time 
intervals and reconstructed object regions, 

Fig. 4 is a flow chart of a filtered back-projection, 

Fig. 5 shows a schematic perspective view of a helical trajectory, a virtual 
20 detector and several beam fans, and 

Fig. 6 shows a further schematic representation of a correlation between 
periods, time intervals and reconstructed object regions. 



25 The computer tomograph shown in Fig. 1 comprises a gantry 1, which is 

capable of rotating about an axis of rotation 14 running parallel to the z-direction of the co- 
ordinate system illustrated in Fig. 1. For that purpose, the gantry 1 is driven by a motor 2 at a 
preferably constant, adjustable angular velocity. A beam source S, for example an X-ray 
tube, is fixed to the gantry 1. The X-ray tube is provided with a collimator arrangement 3, 

30 which from the radiation produced by the beam source S extracts a conical beam bundle 4, 
i.e. a beam bundle that has a finite extent other than zero both in the z-direction and in a 
direction perpendicular thereto (i.e. in a plane perpendicular to the axis of rotation). 
Alternatively, a fan-shaped beam could be used. 
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The beam bundle 4 passes through a cylindrical examination region 13, in 
which a periodically moving object (not illustrated) is located. In this exemplary embodiment 
this object is a beating heart, which performs proper motion and is possibly additionally 
moved back and forth by respiratory motion of the patient. In other embodiments, other 
5 periodically moving body parts, such as the liver, brain or arteries, or periodically moving 
technical objects could alternatively be irradiated. 

After passing through the examination region 13, the beam bundle 4 impinges 
on a detector unit 16 fixed to the gantry 1, the detector unit having a detector surface 
comprising a plurality of detector elements, which in this embodiment are arranged matrix- 

10 form in rows and columns. The detector columns extend preferably parallel to the axis of 
rotation 14. The detector rows are located in planes perpendicular to the axis of rotation, in 
this embodiment on an arc of a circle around the beam source S (focus-centered detector). In 
other embodiments they could alternatively be of a different form, for example, they could 
describe an arc of a circle about the axis of rotation 14 or be linear. Each of the detector 

1 5 elements on which the beam bundle 4 impinges supplies in each position of the beam source 
a measured value for a beam from the beam bundle 4. If in other embodiments a fan-shaped 
beam bundle is used, then the detector unit could alternatively have just a single row of 
detectors. 

The angle of aperture of the beam bundle 4 denoted by a max determines the 
20 diameter of the object cylinder, within which the object to be examined is located during 
acquisition of the measured values. The angle of aperture is here defined as the angle that a 
ray lying in a plane perpendicular to the axis of rotation 14 at the edge of the beam bundle 4 
encloses with a plane defined by the beam source S and the axis of rotation 14. The 
examination region 13 or rather the object, or the patient support table, can be displaced by 
25 means of a motor 5 parallel to the axis of rotation 14 and the z-axis. Alternatively and 

equivalently, the gantry could be displaced in that direction. If the object is a technical object 
and not a patient, the object can be rotated during an examination, whilst the beam source S 
and the detector unit 16 remain stationary. 

By means of the motors 2 and 5, the beam source S and the detector unit 16 
30 are able to describe a trajectory relative to the examination object 13, this trajectory running 
on a notional cylinder surface. This trajectory can run, for example, helically, when both 
motors are in operation. If, however, the motor 5 for advance in the direction of the axis of 
rotation 14 is idle, and the motor 2 allows the gantry to rotate, a circular trajectory is 
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produced for the beam source S and the detector unit 16 relative to the examination region 
13. In the present embodiment, the helical trajectory will be considered. 

During acquisition of the measured values, the cardiac motion is recorded in 
known manner by means of an electrocardiograph 8. For that purpose the thoracic region of a 
5 patient is connected by means of electrodes (not illustrated) to the electrocardiograph 8. 

Alternatively, the pulse could be used as the movement signal describing the cardiac motion. 
In other embodiments, especially in the case of other moving objects, the motion of the 
object can be followed by means of other movement signals. Thus, for example, in the case 
of a technical object that is periodically moved by a motor, a signal of the motor is used as 

1 0 movement signal. 

In the present embodiment, it is assumed that the patient is not breathing 
during the measurement. The respiratory motion can therefore be disregarded. Alternatively, 
the respiratory motion could be measured, for example, using a deformable abdominal belt 
that is connected to a respiratory motion-measuring device. 

15 The measured values acquired from the detector unit 16 are fed to a 

reconstruction unit, especially a reconstruction computer 10, which is connected to the 
detector unit 16, for example, via a contactlessly operating data transmission (not illustrated). 
In addition, the electrocardiogram is transmitted from the electrocardiograph 8 to the 
reconstruction unit 10. The reconstruction unit 10 reconstructs the absorption distribution in 

20 the examination region 13 and reproduces it, for example, on a monitor 1 1. The two motors 2 
and 5, the reconstruction unit 10, the beam source S, the electrocardiograph 8, the 
transmission of the measured values from the detector unit 16 to the reconstruction unit 10 
and the transmission of the electrocardiogram from the electrocardiograph 8 to the 
reconstruction unit 10 are controlled by the control unit 7. 

25 In other embodiments, the acquired measured values and the measured 

electrocardiograms for reconstruction can first be fed to one or more reconstruction 
computers, which forward the reconstructed data, for example, via a fiber optic cable, to an 
image-processing computer. 

Fig. 2 shows the sequence of a measuring and reconstructing method that can 

30 be implemented with the computer tomograph as shown in Fig. 1 . 

After initialization in step 101, the gantry rotates at an angular velocity that in 
this embodiment is constant, but can alternatively vary, for example, in dependence on time 
or on the position of the beam source. 
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In step 102, the examination region or rather the object, or the patient support 
table, is displaced parallel to the axis of rotation and the radiation of the beam source S is 
switched on so that the detector unit 16 is able to detect the radiation from a plurality of 
angular positions. At the same time, or even before the beam source S is switched on, the 
S electrocardiograph 8 is activated, so that an electrocardiogram 21 (Fig. 3) is measured 
simultaneously. 

In the present embodiment, the beam source S moves on a helical trajectory 
relative to the examination region 13. Alternatively, the beams source S could move along an 
arbitrary trajectory running on a notional cylinder surface, for example, along a circular 
10 trajectory. 

In step 103, periods T1...T7 are determined by means of the electrocardiogram 
21 (see Fig. 3). A period in this embodiment is defined by the distance in time between two 
adjacent R-peaks 27 of the electrocardiogram 21. Furthermore, initially in each period 
Ti . . .T7 a time interval Ati . . . At7 is determined. Although the index for the periods and the 

15 time intervals runs only from 1 to 7, the number of periods and time intervals is naturally not 
restricted to this number. On the contrary, the method according to the invention can be 
carried out with an arbitrary number of periods and time intervals. 

The initial determination of the time intervals Ati . . . At 7 can be effected in 
different ways. For example, time intervals of which it is known that the heart moves less in 

20 these time intervals than in other time intervals can be predetermined. Thus, it is known that a 
human heart moves relatively little in a time interval that lies between 65% and 85% of the 
time range between two adjacent R-peaks. 

Alternatively, for initial determination of time intervals in which an object 
moves as little as possible, in the first instance a 4D-data set with a low resolution (e.g. 20 x 

25 20 x 20 cm 3 represented by 64 3 voxels) could be reconstructed, for example by means of a 
filtered back-projection. The 4D-data set contains several CT images of the region of the 
object to be reconstructed (FOV - field of view). Each of these CT images was reconstructed 
exclusively using measured values whose acquisition instants lie in time intervals that are 
identically arranged within the respective period. That is to say, the first CT image has, for 

30 example, been reconstructed exclusively using measured values whose acquisition instants lie 
in time intervals whose midpoints are arranged at a point denoted by 5% RR within the 
respective period. The second CT image has been reconstructed, for example, exclusively 
using measured values whose acquisition instants lie in time intervals whose midpoints are 
arranged at a point denoted with 10% RR within the respective period, and so on. 
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The time intervals that contain acquisition instants of the measured values with 
which the respective CT images have been reconstructed can therefore be arranged staggered 
with respect to one another, for example, by 5% RR, in the respective period. The expression 
"x % RR" in this case denotes a phase point at the location t = tR +0.0 XxA (rr, in which t R is 
5 the instant of the first R-peak of the respective period and AtRR is the distance in time between 
the two R-peaks of the respective period. A phase point is therefore allocated to each CT 
image. Next, CT images having adjacent phase points (e.g. 5% RR and 10% RR) could then 
be compared with one another. For that purpose, for example, absolute differences of 
corresponding voxels of the CT images having adjacent phase points could be formed and 

10 added up, the CT images being more similar, the smaller is the resulting sum. The more 
similar two such compared CT images are, the less the object has moved between the 
corresponding phase points. The time intervals could therefore initially be set in the regions 
of a period in which two CT images having adjacent phase points are as similar as possible. 
If, for example, the CT images having the phase points 70% RR respectively 75% RR are 

15 most similar, then in each period the respective time interval could lie between 70% and 75% 
RR. 

In addition, in step 103, in each initially fixed time intervals a weighting 
function 22 is defined. The weighting function serves to weight each measured value in a 
time interval before it is used for reconstruction of the object. The weighting function 

20 preferably progresses so that the measured values gain increasingly in importance the more 
they are arranged in the middle of the respective time interval. 

The weighting function could alternatively be designed so that measured 
values that lie in a time interval are multiplied by 1 and measured values that lie outside the 
time intervals are multiplied by 0. 

25 In step 104, adjacent periods are considered. At the first execution of step 104, 

first of all the chronologically first time interval Afi of the time intervals initially fixed in step 
103 is determined, which contains acquisition instants of measured values whose 
corresponding beams have irradiated the portion of the object to be reconstructed. The period 
T\ y in which the chronologically first time interval At\ lies, and the subsequent period T 2 form 

30 the first period pair. When step 104 is next executed, the second period pair would be formed 
from the chronologically second period T2 and the chronologically third period 73, and so on. 
In this way, on each execution of step 104, corresponding to the chronological time sequence 
in each case a period pair is formed with periods adjacent in time. 
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The periods of a period pair and the corresponding time intervals are denoted 
in the following description by T\, T 2 and At\, At 2 respectively. These reference numerals, in 
conjunction with Fig. 3, serve to illustrate the following steps and do not, of course, mean 
that the following steps are restricted only to the chronologically first two periods and time 
5 intervals. 

Furthermore, in step 102 the subregion 23 of the FOV traversed both by beams 
whose acquisition instants lie in the time interval A/| of the one period T\ of the period pair 
and by beams whose acquisition instants lie in the time intervals At 2 of the other period T 2 of 
the period pair is determined. Since the acquisition geometry, the acquisition times and the 
initial arrangement of the time intervals are known, this subregion 23 can be determined by 
(simple) geometric considerations or by computer simulations. 

In step 105, the subregion 23 determined in step 104 is reconstructed only with 
measured values whose acquisition instants lie in the time interval At\ of the one period T\ of 
the current period pair, so that a first intermediate image is produced. This reconstruction can 
be carried out using known reconstruction techniques, such as filtered back-projection or 
iterative methods (e.g. ART - algebraic reconstruction technique), wherein, before the 
reconstruction, the measured values whose acquisition instants lie in the time interval A/i are 
weighted in accordance with the weighting function 22. The filtered back-projection is 
explained in more detail below in conjunction with the final reconstruction of the entire CT 
image (see Fig. 4). 

In step 106, a second intermediate image of the subregion 23 determined in 
step 104 is reconstructed using only measured values whose acquisition instants lie in the 
time interval At 2 of the other period T 2 of the current period pair. The reconstruction of this 
intermediate image is also not restricted to specific reconstruction methods. Moreover, here 
too the measured values are weighted in accordance with the weighting function 22 prior to 
reconstruction. 

The reconstruction of the intermediate images can be carried out at a clearly 
lower resolution (e.g. 64 3 voxels for a region of 20 x 20 x 20 cm 3 to be reconstructed) 
compared with the later final reconstruction of the CT image. 

A comparison of the two intermediate images takes place in step 107. For that 
purpose, a similarity measure is applied to the two intermediate images. Many similarity 
measures that compare two images with one another and produce a similarity value 
dependent on the degree of similarity are known. According to the invention, any such 
similarity measure can be used. Known similarity measures use, for example, differences or 
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correlations. The similarity measure a can be formed, for example, by a mean value of the 
absolute differences of corresponding voxels: 

° m jfl£K- v v\ wilhi-1, (1) 

Here, TV is the number of voxels in one of the intermediate images, V\ t \ is the 
image value of the i- A voxel of the first intermediate image and V2 t \ is the image value of the 
i-* voxel of the second intermediate image. Furthermore, the mean quadratic deviation of 
corresponding voxels could also be used as similarity measure: 

^TfZVfW^ with J-l * (2) 

If the similarity measure yields a similarity value that is less than a 
predetermined threshold value, then step 109 is continued. Otherwise step 108 follows. 

In step 108, the time interval A/2 of the later period T2 of the period pair is 
modified. It can be shifted, the width of the time interval can be reduced or enlarged and/or 
the course of the weighting function can be varied. For example, upon each execution of step 
108 the time interval A/2 can be shifted in one direction by a pre-determinable distance of 
time. Thus, at the first execution the time interval A/2 can be shifted by once the distance of 
time in the direction of shorter times, during the next execution by twice the distance of time 
in the direction of longer times, then by three times the distance of time in the direction of 
shorter times and so on. What matters is that upon each execution the time interval A/2, and 
possibly also the weighting function, are modified. 

After the time interval A/2 of the later period Ti of the period pair has been 
modified, in step 106 the subregion 23 determined in step 104 is reconstructed once again 
using measured values whose acquisition instants lie in the modified time interval A/2. The 
resulting new second intermediate image is compared in step 107 by means of a similarity 
measure with the first intermediate image. When the similarity measure yields a similarity 
value that lies above the similarity threshold, then the time interval A/2 of the later period T2 
of the period pair is modified once again in step 108. This time interval A/2 is therefore 
modified until the similarity measure yields a similarity value that is less than the similarity 
threshold. The outcome of this is that, after the similarity threshold is reached, the two 



WO 2005/055828 



11 



PCT/IB2004/052567 



intermediate images are relatively similar. This means that the object states that are 
represented by the two intermediate images are also relatively similar. As will be explained 
in detail below, this leads to a reduction in the movement artifacts in the finally reconstructed 
CT image. 

5 Alternatively, in step 106 the subregion 23 determined in step 104 could also 

be reconstructed repeatedly using measured values from different time intervals At 2 of the 
later period T2 of the period pair. The different time intervals Atz can be arranged close to the 
initial time interval but slightly staggered with respect thereto. The several resulting second 
intermediate images could then each be compared in step 107 with the first intermediate 

0 image by means of the similarity measure, wherein that particular second intermediate image 
(or rather that particular time interval A/2 of the later period Ti) that together with the first 
intermediate image leads to the smallest similarity value, is selected. 

The invention is not restricted to the determination described here of the time 
intervals by means of period pairs. On the contrary, in place of step 104 to 109 it is possible 

5 according to the invention to apply any method that adjusts the time intervals in the periods 
so that a similarity measure applied to several intermediate images of the same subregion is 
minimized, wherein different intermediate images have been constructed using measured 
values from different periods and from only one period in each case. For that purpose, all 
subregions that are reconstructable using measured values from different periods can be 

0 determined. Since an intermediate image is reconstructed for each period, the number of 
intermediate images of a subregion is equal to the number of periods whose time intervals 
contain measured values with which the corresponding subregion is reconstructable. The time 
intervals in the periods are now modified until a similarity measure applied to intermediate 
images of the modified time intervals is minimized. This method can be carried out for each 

5 subregion that is reconstructable using measured values from different periods. 

Step 109 tests whether all periods containing acquisition instants of beams that 
have traversed the FOV have already been taken into account. If that is not the case, then step 
104 is continued. Otherwise, step 110 follows. 

In step 1 10, a CT image of the entire FOV is reconstructed, wherein the 

0 measured values are beforehand weighted in accordance with the respective weighting 
function 22, if this has not already taken place, and wherein only measured values whose 
acquisition instants lie in the time intervals determined in steps 103 to 109 are used. Since the 
time intervals have been determined in such a way that the intermediate images that represent 
the same subregion of the object, yet have been reconstructed with measured values from 
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different periods, are very similar, the object states that are represented by these intermediate 
images are also very similar. If now, as in step 1 10, measured values from different periods 
are used for reconstruction of such a subregion, then the movement artifacts in these 
subregions are reduced in comparison with known methods. 
5 In the reconstruction of periodically moving objects it is naturally 

advantageous if the time resolution in the CT images is as high as possible. A high resolution 
can be achieved by selecting time intervals in the periods that are as narrow as possible. In 
order still to obtain a sufficient number of measured values for the reconstruction of a good 
quality CT image, each subregion of the FOV must be traversed by beams whose acquisition 

10 instants lie in different periods. In the case of known methods, this leads to marked 
movement artifacts, since as a rule the same object state is not found in every period. 
According to the invention, the time intervals in the periods are, however, selected so that the 
above-mentioned intermediate images and hence also the corresponding object states are as 
similar as possible. The method according to the invention therefore leads to reconstructed 

15 CT images having a high time resolution and reduced movement artifacts in comparison with 
the prior art. 

According to the invention, any reconstruction method can be used in step 110 
for the reconstruction of the CT image. For preference, a kind of filtered back-projection is 
carried out, which is illustrated by means of a flow chart in Fig. 4. 

20 For reconstruction, in step 201 the measured values are re-grouped in parallel. 

The measured values are re-sorted and re-interpolated by the parallel re-grouping as though 
they had been measured with a different beam source (an expanded beam source, which is 
arranged on a part of a helix and is able to emit beam fans parallel with one another) and with 
a different detector (a flat, rectangular "virtual detector" containing the axis of rotation 14). 

25 This is explained in detail with reference to Fig. 5. The reference numeral 17 

here denotes the helical trajectory from which the beam source irradiates the examination 
region. A fan-shaped beam bundle 43, the beams of which run in a plane containing the axis 
of rotation 14, is emitted from the beam source position So. The conical beam bundle, which 
is emitted from the beam source in the position So, can be envisaged as comprising a plurality 

30 of flat beam fans that are located in planes parallel to the axis of rotation 14 and intersect at 
the beam source position So- Fig- 5 shows just a single one of these beam fans, namely the 
beam fan 43. 

In addition, Fig. 5 illustrates yet further beam fans 41, 42 and 44, 45, which 
are parallel to the beam fan 43 and lie in planes parallel to one another and to the axis of 
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rotation 14. The associated beam source positions £2, 5-i and S\ 9 S2 are assumed by the beam 
source S before respectively after it has reached the beam source position 5 0 . 

The beam fans 41 to 45 form a group and define a beam bundle 70 having a 
tent-like form. A group of beam fans is called a projection. For each projection a rectangular 
5 virtual detector 160 is now defined which lies in a plane that contains the axis of rotation 14 
and is oriented perpendicular to the parallel beam fans of a projection. The corner points of 
the virtual detector 160 are the penetration points through that plane of the beams that from 
the outer beam source positions meet the opposing portion of the helix. For the ray bundle 70 
in Fig. 5, 5.2, and 52 are the outer beam source positions. Detector elements arranged in 

10 accordance with Cartesian geometry, that is, rows and columns, on which the measured 
values are re-interpolated, are defined on the rectangular detector 1 60. 

In step 202, the measured values associated with the individual beams are 
subsequently multiplied with a weighting factor that corresponds to the cosine of the cone 
angle of the respective beam. The cone angle of abeam is the angle that this beam encloses 

15 with a plane that is oriented perpendicular to the axis of rotation 14. If said angle is small, 
then the cosine of the angle is substantially equal to 1 , so that step 202 can be omitted. 

In step 203, a one-dimensional filtering with a transmission factor increasingly 
linearly with the spatial frequency is applied to the measured values. For that purpose, in 
each case consecutive values in the direction perpendicular to the axis of rotation 14, that is, 

20 along a line of the detector 160, are used. This filtering is performed along each line of the 
virtual detector for all groups of beam fans. 

In another embodiment, the parallel re-grouping could be omitted. Then, as is 
known, since the detector unit is e.g. curved in an arc around the beam source or around the 
axis of rotation, the filtering must be modified. 

25 In step 204, a voxel V(x) inside the FOV is determined. Then, in step 205, a 

projection, that is, a group of beam fans, which has not yet been used for reconstruction of 
the voxel V(x) and whose acquisition instant lies in one of the above-determined time 
intervals, is selected. If no beam of the projection runs centrally through the voxel F(x), then 
it is ascertained at which point a central beam would meet the detector surface. The 

30 associated measured value is then calculated by interpolation of the measured values of 
adjacent beams. The measured value that can be allocated to the beam of the projection 
passing the voxel, or rather the corresponding measured value obtained by interpolation, is 
accumulated in step 206 on the voxel V(x). Step 208 tests whether all projections have been 
taken into account. If this is not the case, then the flow chart branches to step 205. Otherwise, 
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step 208 tests whether all voxels V(x) in the FOV have been passed through. If this not the 
case, then step 204 is continued. If, on the other hand, all voxels V(x) in the FOV have been 
passed through, then the absorption in the FOV as whole is determined and the computer 
tomography method according to the invention ends with step 1 1 1 (see Fig. 2). 
5 In the reconstruction of the CT image, steps 201 to 203 can be omitted if the 

measured values have already been appropriately processed during the reconstruction of the 
intermediate images in steps 105 and 106. 

The method according to the invention is independent of whether the entire 
FOV comprises subregions that are traversed by beams whose acquisition instants lie in 

10 different periods, or of whether only parts of the FOV have these subregions. The latter case 
is illustrated in Fig. 6. 

In Fig. 6 a subregion 24 of the FOV has been reconstructed with measured 
values whose acquisition instants lie in a time interval At'i of the period T\. An adjacent 
subregion 26 of the FOV has been reconstructed with measured values whose acquisition 

15 instants lies in the time interval At'2 in the following period 72. Between these two 
subregions 24, 26 there is a subregion 25, which is traversed both by beams whose 
acquisition instants lie in the time interval At'i and by beams whose acquisition instants lie in 
the time interval At* 2 . In order to minimize movement artifacts in the subregion 25, the time 
intervals At*2, as described above in connection with steps 106 to 108, is modified, so that the 

20 similarity measure applied to intermediate images of the subregion 25, which in each case 
have been reconstructed exclusively with measured values from either the time intervals At*j 
or the time intervals At' 2 , yields a minimized similarity value. The intermediate images are 
hence very similar, which leads to a reduction in movement artifacts in the CT image of the 
entire FOV to be reconstructed. 
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